Frontiers in Systems Biology

Analysis of circadian patterns in time-series data requires computational methods that can accommodate several factors, including variable sampling resolution, replicate number, and missing values. Most existing tools simplify rhythmicity to a strict dichotomy based solely on a single p-value threshold. This leads to a level of uncertainty that affects many biological targets. We developed CODAC (Circadian Oscillation Detection Analysis and Comparison), a framework that integrates nonlinear constrained optimization with a multicriteria rhythmicity classification scheme to evaluate rhythmic patterns without relying on a single statistical cutoff. This approach allows CODAC to identify and exclude medium-confidence rhythms rather than force them into a rhythmic/arrhythmic dichotomy. CODAC comprises four modules: (i) CODAC_single estimates rhythmicity within a single group; (ii) CODAC_flex extends this to identify distinct waveform types within one group; (iii) CODAC_compare performs pairwise comparisons across two or more groups to detect rhythmic or arrhythmic changes; and (iv) CODAC_multi handles more complex designs involving multiple-group comparisons. Using in silico simulations and public transcriptomic datasets, we show that CODAC performs comparably to established methods while providing additional flexibility for rhythm classification and comparison. Taken together, CODAC provides a flexible and open-source package for circadian timeseries analysis with automated visualization tools.

Protein-protein interactions (PPIs) play a crucial role in enabling proteins to carry out their functions within various biological processes (Hui et al., 2003). Since the introduction of the yeast two-hybrid (Y2H) method for PPI detection in 1989 (Fields and Song, 1989), the identification of PPIs has become a significant focus in modern biological research. PPI goes beyond examining individual proteins, allowing researchers to establish a comprehensive network that regulates biological processes. Rice, as a key model organism in plant biological studies, has been at the forefront of PPI research. In 2008, prominent rice scientists in China called for concerted efforts to define a comprehensive protein-protein interaction network experimentally, which aimed to facilitate the prediction of the functional mechanisms operating throughout a plants lifecycle (Zhang et al., 2008). With efforts for 2 decades, the experimentally identified rice PPIs have reached over ten thousand. Several public databases have been established to systematically collate and store PPIs, including STRING (Szklarczyk et al., 2019), BioGRID (Oughtred et al., 2020), IntAct (del Toro et al., 2022), PRIN (Gu et al., 2011), RicePPINet (Liu et al., 2017) and RiceNet v2 (Lee et al., 2015). However, most PPI datasets in rice stem from computational predictions, while experiment-based rice PPI datasets are fragmented due to the lack of systematic profiling at the rice PPIome level, which largely hinders information sharing in the rice research community. To bridge this gap, we constructed the Port of Protein-Protein Interactomes (POPPIN; https://riceome.hzau.edu.cn/poppin/), an integrated database dedicated to sharing experimentally verified PPIs and functional clues in rice. Empowered by high-throughput PPIome profiling technologies and text mining assisted by a large language model (Huang et al., 2025; Liu et al., 2025), POPPIN currently has deposited over 150,451 pieces of rice PPI-related information. Additionally, POPPIN provides detailed protein information, including GO annotations, subcellular localizations, domains, trait ontology (TO) information, and hyperlinks to external biological databases. Through offering a user-friendly web interface for search and dynamic network visualization, POPPIN serves as the first large-scale, experiment-based database for searchable PPIs in rice, and has the potential to be extended to other species under this structural framework.

Selecting stochastic gene-expression models from single-cell counts requires accurate parameter inference and efficient model selection. Likelihood methods in count space can be costly when full stationary count distributions are unavailable, whereas approximate methods may lose accuracy. Probability generating functions (PGFs) offer a compact analytical alternative, but existing PGF workflows are generally not likelihood based and therefore rely on computationally intensive cross-validation. We develop a likelihood-based PGF framework for both tasks. Correlated empirical PGF values are used to construct a Gaussian quasi-likelihood for parameter inference and PGF-based Bayesian information criterion (BIC) for model selection. We show that the empirical PGF is exactly unbiased and that the parameter estimator is consistent, converges at the inverse-square-root sample-size rate, and is first-order asymptotically unbiased. For large samples and a uniquely preferred model, PGF-BIC selects the same model as leave-one-out cross-validation in PGF space.

Spatial transcriptomics keeps measurement of gene expression while preserving spatial context, yet traditional analysis methods face challenges in computational efficiency, biological interpretability, and autonomous discovery. This project presents a framework solving these issues through three parts: (1) an ensemble clustering system achieving 66.7% improvement over baseline average and 23.9% over best single method with silhouette score of 0.540 and statistical significance (p = 0.0032, Cohens d = 1.82); (2) a knowledge-based clustering framework that annotates 88.6% of cells across 8 ovarian cell types using 428 marker genes; and (3) a GPT-4o-mini-powered autonomous agent that generated 3 biological hypotheses with validations.

Long-term in vivo transcriptomic time courses are costly, limiting assessment of chronic molecular responses from short studies. We developed a pharmacokinetics-informed transcriptomic ordinary differential equation model (PKT-ODE) that links an oral pharmacokinetic profile and Hill drug-effect function to first-order turnover of co-expression modules. The model was fitted to rat liver responses to fenofibrate at three doses in Open TG-GATEs through day 8. At the held-out day-29 endpoint, PKT-ODE achieved Pearson r = 0.960 and mean squared error (MSE) = 0.148. In this dataset, these values achieved lower prediction error and higher correlation than four statistical baselines and validation-selected linear and multilayer-perceptron transition models. Literature-curated peroxisome proliferator-activated receptor target genes occurred only in modules with positive fitted drug effects. These results provide a proof of concept for pharmacokinetics-informed transcriptomic extrapolation; cross-compound, cross-organ and alternative-regimen performance remain to be tested.

Metastasis involves context-dependent molecular interactions in which non-coding RNAs, particularly miRNAs and circRNAs, play important regulatory roles. However, existing computational approaches generally do not jointly represent cancer type, metastatic event, and cancer-specific metastatic context. We developed a context-aware multi-task heterogeneous graph neural network (GNN) for predicting ncRNA associations with cancer types and metastatic events. The framework integrates multiple biological repositories into a heterogeneous graph representing ncRNAs, cancers, metastatic event types (METs), and cancer-specific metastatic instances (CSMIs). The model performs six link-prediction tasks using a hierarchical transformer-based encoder and multi-relational TuckER decoder. Across ten independently initialized runs evaluated on the RNA-group-disjoint held-out test set, the model achieved a global AUROC of 0.8801 {+/-} 0.0118 and an F1 score of 0.8260 {+/-} 0.0071. All three ablation variants yielded lower AUROC, with the largest reduction under independent task training. Case studies in pancreatic cancer, colorectal cancer, and hepatocellular carcinoma provided disease-level, event-level, and expression-based support, respectively, for top-ranked candidate associations. The framework enables context-specific prioritization of ncRNA-cancer-metastasis associations for experimental evaluation.

Light history leaves persistent changes in circadian period, but where this history is stored remains unresolved. Suprachiasmatic nucleus (SCN) network models have often approached photoperiodic encoding through phase organization or coupling strength. We computationally tested slow adaptation of subregion-specific intrinsic periods as an alternative memory mechanism. The model asymmetrically couples dorsal (D) and ventral (V) SCN oscillators and adds a systemic oscillator (X) representing putative circadian feedback present in vivo but lost ex vivo. With a single parameter set, period adaptation captured the direction and approximate magnitude of behavioral aftereffects across photoperiod and T-cycle conditions. Adapting coupling strength instead of period failed to reproduce the V-leading-D phase order reported after T22. Removing systemic feedback preserved the photoperiod-dependent period ordering but inverted the T22 and T26 aftereffects, an inversion that matched SCN explant observations. The model also yielded distinct D-V phase organization for each of 18:6 LD, T23, and T25. These results suggest that subregion-specific period plasticity provides a parsimonious substrate for encoding light history, while the dependence on systemic feedback indicates that behavioral period may not be a readout of the SCN alone. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=160 SRC="FIGDIR/small/743784v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@d48266org.highwire.dtl.DTLVardef@1bd15bdorg.highwire.dtl.DTLVardef@de3d9forg.highwire.dtl.DTLVardef@9fbc39_HPS_FORMAT_FIGEXP M_FIG C_FIG A model with dorsal period adaptation and phenomenological systemic feedback accounts for behavioral and explanted SCN aftereffects. (A) During T22 entrainment, dorsal (D), ventral (V), and systemic (X) oscillators remain phase-locked. After release into constant darkness, systemic coupling maintains a unified in vivo rhythm, whereas removing X feedback in the explant simulation allows the D-V network to express a distinct ex vivo period aftereffect. (B) The SCN model is modeled as an asymmetrically coupled attractive-repulsive oscillator network with stronger photic input to V. Light history is encoded via plasticity of the intrinsic period in D, while X represents putative systemic circadian feedback available in vivo and lacking direct photic input. (C) The model reproduces concordant period changes in behavior and SCN explants across photoperiods, but opposing period changes following T-cycle entrainment.

Cross-dataset generalization of EEG-based classification under weak, proxy-derived labels remains an open problem for altered-states research. We present a reproducible eight-dataset alignment pipeline that maps eight heterogeneous EEG corpora (712,832 windows; 697,906 with valid labels) to a common 14-channel EPOC+ montage with 63-dimensional spectral features, and we recover the real 1-9 arousal self-assessments for MAHNOB-HCI from session.xml metadata. As a benchmark, Random Forest classifiers are trained on seven source domains and evaluated on the held-out target under both zero-shot and 20%-participant few-shot calibration. The benchmark exposes two concrete methodological pitfalls rather than a performance result: (i) per-class recall shows every target collapsing to a single majority class, and (ii) a within-dataset upper-bound experiment (Table 3) shows that six of eight proxy label sets sit at or below three-class chance even when trained and tested on the same dataset, so the cross-dataset failure is a label-validity problem rather than a transfer-method problem. Across the eight targets (20 seeds, 8,000 evaluation windows per target), zero-shot accuracy averages 36.85% (95% CI 34.40-39.30) and calibrated 43.76% (41.77-45.75), but balanced accuracy stays at 33.01-35.62% (Cohen's kappa <= 0.068), i.e. at chance. The +6.91pp mean change is driven almost entirely by a single target, ds006437 (6.31% -> 60.60%): the median paired change across all 160 seed-pairs is 0.00pp, and after Holm-Bonferroni correction only ds006437 and ds004572 remain significant, the latter with a practically null effect (+0.39pp). Balanced accuracy stays between 33.01% and 35.62% and Cohen's kappa at 0.009 +/- 0.032, i.e. at or barely above three-class chance, while per-class recall shows six of eight targets collapsing to Deep (96.7-100% recall) and two to Light (68.5-99.1%). The collapse persists under SMOTE oversampling, under an EEGNet-v4 deep-learning baseline, and under CORAL and AdaBN feature alignment, which locates the bottleneck in proxy-label validity and class overlap in the feature space rather than in classifier capacity. We position this work as a preliminary methodological study: its contribution is a reproducible eight-dataset alignment pipeline, recovered MAHNOB-HCI arousal self-assessments, a quantitative estimate of split-leakage inflation, and a transparently reported negative result rather than a performance claim.

Subject-independent affect regression from physiological signals remains challenging because emotional responses vary substantially across individuals, widely used datasets provide only coarse trial-level annotations, and heterogeneous physiological modalities may not contribute reliably when treated as if they were interchangeable predictors. We have developed AffectRoute, a protocol-conditioned subject-independent affect regression that is conditioned on the protocol and assigns separate predictive functions to information from the population, electroencephalography (EEG), and peripheral physiological signals (PPS). First, a source-population prior establishes a trial-level affective anchor using only the data from source participants. TrajBridge then combines an EEG representation that is supervised by REFED for participant-specific adjustments with temporal structure obtained from the continuous REFED annotations in order to create a weakly supervised segment-resolved pseudo-trajectory and to establish a frozen trial-level baseline. PhysioRoute next reduces the remaining error by breaking down the residual correction into a source-derived direction, which is estimated from the out-of-fold residuals within the source group, and a channel-specific magnitude derived from the PPS. When evaluated on DEAP and DREAMER using a leave-one-subject-out approach at the participant level, AffectRoute showed consistent step-by-step improvements in both the mean absolute error and the concordance correlation coefficient. A method that relied solely on the source data was clearly worse than PhysioRoute, showing that the final improvement cannot be accounted for by transferable source residual regularity alone. Conventional alternatives to fusing the PPS were also found to be consistently less effective, although analyses at the channel level and with a leave-one-channel-out design showed that the peripheral contributions are axis-dependent yet distributed across channels. These results indicate that structured residual inference is an effective alternative to unrestricted multimodal fusion for subject-independent affect regression.

Recent advances in AI co-scientists have brought LLM agents into closed-loop experimental design. However, whether these agents use feedback from earlier rounds to revise subsequent experimental decisions remains unclear. We address this question with PerturbTrace, which evaluates each round-to-round transition through Feedback-to-State, State-to-Action, and Action-to-Outcome. These stages assess whether feedback is reflected in the agent's rationale and perturbation-selection strategy, whether the stated strategy guides the next perturbation batch, and whether that batch yields more hits than expected under random sampling. We evaluate four LLM agents on 17 screen-derived tasks and compare them with random selection, active learning, and LLM-guided Bayesian optimization baselines. Each agent outperforms the strongest non-agent method on at least 15 of the 17 tasks, yet controlled evaluations across six tasks show no consistent advantage from true feedback over random or no feedback. Among 576 transitions under true or random feedback, only 43 (7.5%) complete the full Feedback-State-Action-Outcome sequence, including 25 under random feedback. These findings show that high final recall does not necessarily indicate effective feedback use. They also highlight the need to evaluate closed-loop scientific agents by both their discovery performance and whether feedback changes their subsequent decisions.

Isotopic Metabolic Flux Analysis (I-MFA) is a standard approach for estimating intracellular metabolic fluxes. I-MFA infers fluxes by comparing simulated and measured metabolite isotopologue distributions (MIDs) of metabolites from isotope labeling experiments. MIDs represent fractional abundances that strictly sum to one for any given metabolite, thus they are inherently compositional data. However, state-of-the-art estimation approaches rely on calculating standard Euclidean distances between MIDs in a non-compositional paradigm, introducing a systemic bias. To resolve this, our study proposes compositional I-MFA. We demonstrate how to construct a meaningful orthonormal basis for MIDs via ordered sequential binary partitioning, which can be used to perform isometric log-ratio (ILR) transformation. As a minimal change to existing I-MFA workflows, we suggest estimating fluxes by minimizing Euclidean distances between ILR-transformed MIDs. We validated this framework against traditional methods using both a toy model and a biologically realistic model, evaluating point estimates, sensitivity across varied true fluxes, and confidence intervals. In the two examples, compositional I-MFA consistently outperformed traditional approaches, reducing mean squared error of flux point estimates by an average of 42.6% and substantially narrowing confidence intervals. We conclude that compositional data analysis significantly improves I-MFA and can be implemented as a simple drop-in replacement for current pipelines. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=156 SRC="FIGDIR/small/742769v1_ufig1.gif" ALT="Figure 1"> View larger version (34K): org.highwire.dtl.DTLVardef@1a53aa4org.highwire.dtl.DTLVardef@ad225aorg.highwire.dtl.DTLVardef@aa430eorg.highwire.dtl.DTLVardef@1880ca_HPS_FORMAT_FIGEXP M_FIG C_FIG HighlightsO_LINew compositional data approach improves metabolic flux estimation. C_LIO_LIThis data transformation requires minimal changes to existing workflows. C_LIO_LIThe new method reduced MSE of flux estimates by 42.6% in two examples tested. C_LIO_LIThe confidence intervals of the estimated fluxes were substantially narrowed. C_LIO_LIEstimation accuracy remained robust across a wide range of metabolic fluxes. C_LI

Mechanistic drug-development programmes often have more biological evidence than they can safely quantify. We developed evidence-constrained mechanistic synthesis (ECMS), a framework that classifies what information each finding contains and converts only that information into restrictions on a family of mechanistic hypotheses. Evidence shifts the frequency of supported events in a reproducible ensemble rather than being converted into unsupported coefficients or probabilities of biological truth. In a chronic spontaneous urticaria (CSU) implementation, a representative, non-exhaustive corpus of 114 atomic findings from 53 sources and 13 public data resources compiled 18 relation/context constraints and a frozen 4,096-hypothesis ensemble. Regimen evaluation was formulated as continuous multi-node target matching: researchers specify desired changes and importance coefficients for modeled nodes, while package-declared controls vary continuously. A deterministic Sobol-to-block-refinement search, validated on all 4,096 hypotheses, reduced target-matching loss by 27.3% relative to the best of 44 deterministic anchors under a prespecified heuristic demonstration profile; changing the objective profile changed the selected control vector without changing the evidence ensemble. A complementary D-only reference analysis localized decision-relevant uncertainty around the mast-cell-to-disease relation, illustrating that mechanistic prioritization depends on the declared objective. ECMS is intended for the pre-calibration stage of drug development: it makes heterogeneous literature computable while keeping evidence, uncertainty and decision preferences distinct.

Control of gene expression by transcription factors (TFs) is a critical mechanism for cells to maintain homeostasis in response to environmental signals. Gene network models that predict regulatory interactions between transcription factors and the genes they control aid in understanding these complex processes. These models are useful as they provide testable hypotheses of regulatory interactions, transcription factor function, and accelerate the study of uncharacterized transcription factors. However, inference of these models is computationally challenging due to the vast quantity of data required given the many possible states of the regulatory network. Microbial genomes encode hundreds of transcription factors, with numerous interactions that require substantial functional genomics datasets to infer. This problem is accentuated in understudied organisms, species that would greatly benefit from an inferred network for biological discovery, where the lack of available data is particularly constraining for effective inference. To address this problem, we have developed GRN-BMuSeR (Gene Regulatory Networks from Bayesian MUlti-SpEcies Regression), a novel multitask approach to gene regulatory network inference that leverages gene orthology between closely related species to improve inference performance. We evaluate its performance on a dataset from the well-studied bacterial species Bacillus subtilis, demonstrating improved performance in multitask settings. Applying the model to simulated data reveals utility in multi-species contexts. Finally, we apply our models to infer GRNs and explore predictions for two hypersaline-adapted archaeal species. We leverage a rich dataset from Halobacterium salinarum to inform the inference of the gene regulatory network of Haloferax volcanii, for which a more limited genomics dataset was available. We generate a large compendium of gene expression data for Hfx.volcanii for GRN inference input. Through exploration of resultant network predictions, we show concordance with known TF functions and discover hundreds of novel TF functional predictions. Moving forward, our results provide a framework to generate testable hypotheses that will serve to guide experimental work and accelerate discovery in these understudied species.

Tumors consist of heterogeneous phenotypic cells, such as normoxic cells, which are highly proliferative, and hypoxic cells, which are less proliferative. Their phenotypic switching depends on tumor microenvironmental factors, such as oxygen and nutrient concentrations supplied by local blood vessels. However, during ongoing angiogenesis, the process of sprouting new blood vessels at the tumor site from pre-existing blood vessels, and how this phenotypic switching affects and impacts tumor growth, remains poorly understood. In this article, we formulate a mathematical model to elucidate the crosstalk between vasculature and tumor cellular heterogeneity during tumor progression. The model results show a strong agreement with the experimental data. Our simulation results demonstrate that ongoing angiogenesis increases tumor growth rate. In addition, we observe that the influence of hypoxic cells on phenotypic switching from normoxic to hypoxic is more pronounced than their influence on the transition from hypoxic to normoxic. Furthermore, we perform a global sensitivity analysis using the Sobol's method to assess the importance of the model's parameters. It highlights that the volume at which blood vessels attain half-maximal rate has the maximum effect on the model.

We introduce LifeSciBench, a benchmark of 750 expert-authored tasks designed to evaluate whether language models can handle realistic life science research work. The majority of existing life sciences benchmarks have a narrow scope or are purely knowledge-based, and therefore fail to capture the complexity of real-world research, which often involves ambiguities and requires the accurate execution of multiple dependent judgment calls. Additionally, almost all existing benchmarks span at best a small collection of subdomains within the life sciences; there is at present no existing life sciences benchmark with both the requisite breadth and depth required to convincingly measure proficiency in real-world professional research settings. LifeSciBench addresses this gap by spanning seven representative scientific workflows and seven life science domains, with each constituent task paired with a human expert-written rubric. Across five frontier and domain-specialized models, GPT-Rosalind performs best, with a task-weighted mean normalized rubric score of 0.576 and a task-weighted response pass rate of 36.1% (response-level values are first averaged within each task, and the resulting task-level values are then averaged with equal weight). LifeSciBench remains unsaturated, with 171 tasks (22.8%) having no observed passing response from any evaluated model and 261 tasks (34.8%) having a best-model pass rate below 20%. LifeSciBench therefore serves as a high-resolution evaluation of practical scientific reasoning and operational decision-making in the life sciences.

Summary: Agentic AI platforms enable the engineering of autonomous workflows but are not designed for experimentation and hypothesis testing. ASAREE (Analytical Sandbox for Agentic AI Research, Engineering, and Experimentation), is an open-source platform to address this gap. ASAREE creates agents, connects to MCP servers and tools, and designs factorial experiments through a visual interface or Python SDK. It records a full provenance trace for every run and routes all model calls through a provider-agnostic bridge that supports local deployments, ensuring data privacy. As a use-case, we use ASAREE to evaluate key design choices in a mutli-agent machine learning pipeline. Across a 2 x 2 x 2 factorial design, more advanced models, greater reasoning effort, and critic agent use significantly increased compute time, token use, cost, and feature count without improving predictive performance. The lowest-cost baseline, Claude Sonnet 5 with medium effort and no critic, achieved the highest mean PR AUC while Claude Opus 5 with extra high effort and a critic agent cost 15.5x more (USD) and ran 13.1x longer while performing worse on average. These findings highlight ASAREE as a robust framework for evaluating agentic system performance and resource efficiency.

Machine learning (ML) models for molecular property prediction are increasingly deployed in drug discovery, yet their adoption in real-world scenarios requires an understanding of the conditions in which a model succeeds or fails. While standardized benchmarks are powerful instruments to measure and unlock progress in ML research, they should not be blindly treated as the end goal. Especially static and retrospective benchmarks, in which no true unknown test set is employed, limit our ability to robustly validate a model's performance. Building on the collective expertise of a cross-industry consortium, we present a model validation framework consisting of five recommendations that would enable the community to move beyond aggregate metrics toward understanding where and why molecular property prediction models fail. We connect evaluation choices to real-world applications and case studies encountered in pharmaceutical research. The framework proposes splitting strategies that mimic realistic distribution shifts and expose common failure modes. We apply the recommended framework to a recently released dataset of absorption, distribution, metabolism, and excretion (ADME) properties. Across two complementary model algorithms, our case studies reveal four distinct failure modes (extrapolation, interpolation, representation, and evaluation), showing that model errors arise not only from distribution shift but also from limitations in molecular representations. Our results show that commonly used evaluation protocols can significantly overestimate performance and may not detect important model failure modes. All software and data are released via https://github.com/srijitseal/polaris.

Drosophila connectomic datasets provide increasingly comprehensive maps of neuronal morphology and synaptic connectivity, offering an unprecedented opportunity to explore the structural organization of its neural circuits. This calls for designing automated tools to interact with connectomic datasets at scale for efficiently exploring structural features embedded in the vast amount of data. Yet the central challenge remains the understanding of the functional logic of neural circuits. In order to understand how elements of the functional logic may emerge from this structural organization, it is critical to (i) characterize the objects in the natural environment in which brain circuits operate, and (ii) formulate how brain circuits represent and process the defined objects in the natural environment. To develop and demonstrate a methodology for these requirements, we focus on the Drosophila looming-evoked escape pathway. We modeled the trajectory of looming objects that are on a collision course (direct-hits) or pass-by the fly (near-misses): their projected images on the retina can be characterized by the solid angle (angular size) and elevation. We then analyzed the pathway's morphology across the OpticLobe, Hemibrain, and FlyWire connectome datasets. By abstracting their sub-neuronal structure and retinotopic organization, we constructed an executable circuit model that maps each structural element to a processing block. We demonstrate that this model separates direct hits from near misses well before the angular size could tell them apart. To accelerate the connectomic analysis step, we developed a Python toolset with an agentic, code-free workspace interface called NeuroGraphBench (NGB). NGB provides four composable morphology-analysis primitives and an AI agent that composes them to interactively respond to natural-language queries aided by visualization on an interactive 3D canvas. Thus, NGB automates tedious and repetitive tasks to enable faster and scalable connectomic exploration, keeping human reasoning, instead of writing code, at the center of an open-ended research inquiry.

Light simultaneously provides phototrophic organisms with energy and with information about environmental time. These two functions need not impose the same response to fluctuations in irradiance: photosynthetic outputs should remain amplitude-sensitive, whereas circadian phase should reject changes that do not alter dawn, dusk, or photoperiod. We formulate this problem for two spectral inputs by decomposing their logarithmic intensities into a common-irradiance coordinate a and a spectral-contrast coordinate r. The contribution of channel i to phase is Qi = ZiGi, where the non-negative gate Gi determines when the pathway is active and the signed phase-response projection Zi determines whether this activity advances or delays the oscillator. For a locked oscillator, robustness to common irradiance together with retained contrast sensitivity requires two non-zero cycle-averaged contributions of opposite sign, A1 [~=] -A2 = 0. Energetic responsiveness is preserved only when the physiological projection of the same inputs is not proportional to their phase projection. A canonical repressilator provides an explicit nonlinear realization of these conditions. Positive gates placed on opposite lobes of its infinitesimal phase-response curve strongly attenuate common-mode phase shifts while preserving contrast sensitivity. A minimal photosynthetic-capacity model then shows how this organization protects temporal alignment under day-to-day irradiance fluctuations. At the largest variability tested, differential routing reduced the mean phase displacement by more than one half and the associated alignment loss by approximately 82%, whereas the resulting production advantage remained small, approximately 0.1%. Thus, multichannel light sensing can stabilize circadian timing without suppressing the energetic response to irradiance. HighlightsO_LIAnalytical routing conditions separate common irradiance from spectral contrast. C_LIO_LIPositive temporal gates can generate opposite signed phase contributions. C_LIO_LIPhase robustness requires a projection distinct from the energetic projection. C_LIO_LIA canonical oscillator provides a constructive illustration of the mechanism. C_LIO_LIThe functional benefit is improved temporal alignment rather than a large growth gain. C_LI

Determining the amount of information a neuron receives per unit of time is key to understanding brain connectivity and how neural networks encode and transmit information. In particular, estimating this information flow by distinguishing between excitatory and inhibitory synaptic contributions is critical to understanding neural network function, as maintaining the excitation-inhibition (E/I) balance regulates neuronal excitability and circuit stability, whereas its disruption can lead to a plethora of brain disorders, including neurodegenerative and psychiatric conditions. However, because synaptic conductances cannot be measured directly, inverse methods are required to infer them from the membrane potential --a readily measurable quantity. Although partial solutions have been proposed, accurately estimating these conductances remains a significant challenge due to the complexity and diversity of the inputs. This is particularly true in the spiking regime, where neurons actively fire. In this work, we introduce a novel computational strategy that combines two critical metrics extracted from the time course of the membrane potential recording: the amplitude of the spike and the interspike interval. By using these quantities, the proposed method enables the accurate separation of excitatory and inhibitory contributions, yielding highly favorable results in the spiking regime. Author summaryQuantifying the continuous stream of inputs a neuron receives is key to understanding brain connectivity. Inside the brain, individual cells must maintain a tight balance between excitation and inhibition (E/I) to process information correctly, as any disruption in this equilibrium can impair its functionality. However, directly measuring the underlying excitatory and inhibitory synaptic conductances is technically challenging, and existing mathematical tools often fail when neurons enter their active firing regime. In this work, we introduce a novel computational strategy designed to extract and separate these time-varying conductances directly from the neurons spiking activity. By dynamically tracking just two accessible metrics - the amplitude of the spikes and the time intervals between them - our algorithm estimates both conductance profiles with high precision. Furthermore, we demonstrate that this procedure is highly robust against realistic experimental noise and data variability, providing an accessible framework that does not require complex hardware or an unfeasible number of repetitive experimental trials. By tracking changes in the E/I ratio of the synaptic input, this method provides an efficient approach to detecting pathological imbalances and understanding how local connectivity shapes cellular functionality.